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APPLIED  TECHNOLOGY  LABORATORY  POSITION  STATEMENT 


This  report  provides  the  documentation  for  a  program  which  demonstrated  a  method  for 
balancing  a  full-size  supercritical  helicopter  tail  rotor  drive  shaft  using  a  quasi-static  data 
acquisition  technique.  The  process  has  been  shown  to  place  the  center  of  gravity  at  all 
stations  within  0.001  inch  of  the  center  of  rotation.  The  balancing  technique  is  available 
for  implementation  in  existing  and  future  tail  rotor  drive  shafts. 
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and  Mr.  Fred  Reed  of  the  Directorate  for  Systems  Engineering  and  Development,  US  Army 
Aviation  and  Research  Development  Command  served  as  project  engineers  for  this  effort. 
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SUMMARY 


This  report  presents  the  results  of  a  balancing  technique  for  flexible  rotors. 
The  technique  was  demonstrated  using  the  YAH-64  Phase  I  Helicopter  Tail 
Rotor  Drive  Shaft  as  the  flexible  rotor. 

The  technique  measures  the  mass  distribution  of  the  rotor  about  its  center  of 
rotation  at  a  series  of  stations  along  its  length.  It  uses  ultrasonic  gaging 
equipment,  a  digital  computer  to  process  the  50,000  data  points  measured  by 
the  equipment,  and  printout  to  define  the  localized  balance  needed  to  bring  the 
local  center  of  gravity  back  to  the  center  of  rotation.  The  process  has  been 
shown  to  place  the  center  of  gravity  at  all  stations  within  0.  001  inch  of  the 
center  of  rotation. 
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PREFACE 


This  report  was  prepared  by  Hughes  Helicopters,  Division  of  Summa 
Corporation,  under  Contract  No.  DAAK5  1  - 78- C-002  1 ,  funded  by  Applied 
Technology  Laboratory,  USAA VRADCOM,  Fort  Eustis,  Virginia.  It  covers 
the  work  performed  during  the  period  September  1978  to  February  1980  and 
is  the  final  technical  report  summarizing  the  activity.  The  contract  statement 
of  work  breaks  down  the  activities  into  six  tasks. 

TASK  1 

A  search  was  made  for  manufacturers  of  ultrasonic  gaging  equipment.  Three 
companies  were  contacted  for  literature.  All  responded,  but  only  one 
(N.D.  T.  Instruments)  was  able  to  delivv.  •  hardware.  It  was  concluded  that 
the  ability  to  obtain  the  measuring  requirements  for  this  program  using 
ultrasonics  was  "state  of  the  art". 

TASK  II 

Design  and  fabricate  the  equipment  required  to  modify  an  existing  Hughes 
Helicopters'  balancing  fixture  (used  during  the  Phase  I  development  of  the 
YAH-64  Helicopter),  for  compatibility  with  the  automatic  data  acquisition  sys¬ 
tem  and  interface  with  the  computer.  Create  the  computer  programs  required 
to  perform  the  calibration,  acquisition,  and  data  functions  analysis.  The  final 
product  is  a  plot  and  tabulation  of  the  necessary  balance  corrections. 

TASK  III 

Select  two  shafts  that  were  manufactured  during  the  Phase.  I  development  of 
the  YAH-64  Helicopter  Program.  The  two  shafts  (P/N  7-113500003)  were 
S/N-ll1,  used  in  the  Phase  I  YAH-64  ground  test  vehicle,  and  S/N-5U,  a  spare 
for  air  vehicles  02  and  03. 

TASK  IV 

Trial  balance  S/N-1U  for  checkout  and  evaluation.  , 

Balance  S/N-5U  witness  by  the  Contracting  Officer's  representative. 

TASK  V 

Conduct  an  economic  analysis  of  this  method  for  balancing  1,000  shafts. 

TASK  VI 

Present  a  Government/ industry  briefing  at  Hughes  Helicopters'  facility, 

Culver  City,  California. 
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INTRODUCTION 


As  designers  become  more  weight  and  vibration  conscious  for  new  aircraft 
designs,  the  introduction  of  power-transmitting  shafts  that  run  at  super¬ 
critical  speed  becomes  increasingly  attractive.  Balancing  such  high-speed 
shafts  so  that  they  run  smoothly  has  long  been  a  problem,  with  the  process 
being  both  difficult  and  costly.  It  traditionally  requires  many  trial  runs  over 
the  rpm  spectrum  to  achieve  good  balance. 

The  technique  described  in  this  report  was  developed  by  Hughes  Helicopters 
(HH)  and  has  been  granted  U.  S.  Patent  No.  4,  170,  896  dated  October  16, 

1979.  The  balance  method  requires  no  trial  runs  to  develop  baseline  data. 

All  required  data  is  obtained  during  a  "one  shot"  pass  along  the  slowly- 
rotating  shaft,  is  automatically  processed,  and  the  operator  is  given  a  graphic- 
display  of  how  much  balance  is  to  be  applied  along  the  shaft.  The  system  is 
shown  schematically  in  Figure  1. 

Figure  1  is  a  schematic  representation  of  the  balancing  process.  The 
ultrasonic  transducer  and  sensors  ride  along  a  track  that  is  parallel  to  the 
shaft's  centerline  of  rotation  and  measure  two  quantities  at  each  of  <150  azi¬ 
muthal  locations  around  the  shaft  for  each  lengthwise  station  along  the  shaft. 
The  "slave"  sensor  measures  the  distance  from  the  transducer  to  the  outside 
surface  of  the  shaft  while  the  "master"  measures  the  thickness  of  the  wall  of 
the  shaft.  The  computer  integrates  this  information  around  the  periphery  of 
the  shaft,  calculates  the  local  center  of  gravity,  calculates  the  amount 
and  location  of  balancing  tape  needed  to  bring  the  local  center  of  gravity  and 
center  of  rotation  into  coincidence,  and  displays  the  balance  requirements  in 
both  analog  and  digital  form.  The  operator  then  applies  the  adhes  ive -backed 
aluminum  tape  to  the  surface  of  the  shaft. 

This  method  is  valuable  not  only  as  a  shaft-balancing  procedure  but  also  as 
a  low-cost  receiving  inspection  tool  to  qualify  the  basic  tube  for  the  important 
parameters  which  make  an  acceptable  supercritical  drive  shaft,  i.  e.  , 
straightness,  wall  thickness  variations,  and  dents.  By  rejecting  tubes  that 
require  large  or  impractical  balance  adjustments,  all  the  labor  of  shaft  manu¬ 
facture  is  saved,  thereby  minimizing  the  rejection  rates  for  a  finished  shaft. 
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INSTRUMENTATION 


The  instrumentation  required  to  conduct  this  program  can  be  categorized 
into  five  basic  areas: 

•  Ultrasonic  gaging 

•  Preamplification 

•  Analog-to-digital  conversion 

•  Shaft  azimuth  encoding 

•  Digital  computation  and  output 
ULTRASONIC  GAGING 


The  ultrasonic  equipment  was  purchased  from  N.D.T.  Instruments,  Inc., 
Huntington  Beach,  California.  It  consists  of  two  specially  prepared 
Nova  Scope  2000  instruments  and  one  ultrasonic  transducer  with  a  water 
coupling  adapter  and  a  mechanical  manipulator.  One  of  the  instruments  acts 
as  a  slave,  while  the  other  performs  as  the  master. 

The  master  unit  provides  the  excitation  pulse  to  the  transduce  r,  and  both  Nova 
Scopes  read  back  the  ultrasonic  signal  response.  The  slave  unit  is  calibrated 
to  provide  data  on  the  "water  path"  or  distance  from  the  transducer  to  the  top 
surface  of  the  specimen  (drive  shaft  surface).  The  master  is  set  up  and 
calibrated  to  provide  the  wall  thickness  data. 

Both  pieces  of  data  are  extracted  simultaneously  from  the  same  ultrasonic 
pulse.  Figure  3  shows  a  typical  ultrasonic  pulse  and  the  interpretation  of 
the  peaks.  The  slave  signal  circuitry  employs  a  blocking  gate  to  delay  the 
timing  circuit.  This  gate  is  set  up  by  the  operator  to  some  short  time  before 
the  first  echo.  The  purpose  is  to  avoid  false  readings  due  to  environmental 
noise.  The  time  between  the  main  pulse  and  the  first  return  echo  is  equated 
to  twice  the  distance  to  the  first  surface  being  measured.  The  instrument 
is  calibrated  by  making  precise  known  changes  in  the  distance  to  the  tirst 
surface,  and  adjusting  the  gain  and  zeroing  potentiometers  on  the  front  ot 
the  Nova  Scope.  The  digital  display  on  the  front  of  the  Nova  Scope  provides 
a  convenient  indication  for  calibration  and  can  be  used  at  any  time  to  indi¬ 
cate  the  parameter  being  gaged. 


The  master  unit  circuitry  is  set  up  to  block  out  the  main  pul_e  and  measure 
the  time  from  the  first  return  echo  to  the  second  return  echo.  This  time  is 
equated  to  the  wall  thickness.  This  unit  is  calibrated  by  using  several  known 
thicknesses  of  aluminum  and  adjusting  the  gain  and  zero  potentiometers  until 
the  digital  display  indicates  the  thickness  of  the  sample  being  gaged.  An 
electrical  connection  to  an  analog  voltage  that  is  proportional  to  the  digital 
signal  displayed  on  the  front  is  provided  at  the  rear  of  each  Nova  Scope  2000. 
This  analog  voltage  is  the  input  to  the  computer. 

PREAMPLIFIERS 

Standard  DC  laboratory  instrumentation  amplifiers  serve  as  preamplifiers 
between  the  analog  output  of  the  Nova  Scopes  and  the  analog -to-dig ital  (A/D) 
converters  of  the  computer  to  produce  a  digital  resolution  at  the  computer 
compatible  with  the  desired  sensitivities;  i.  e.,  less  than  0.001  inch  for 
eccentricity  measurement  and  0.0001  inch  for  the  wall  thickness  measurement. 

DIGITAL  COMPUTER  A/D 

A  Hughes  Helicopters'  laboratory  computer  is  used  for  this  program.  It  is  a 
Hewlett  Packard  1000  computer  with  a  HP  2313  Data  Acquisition  System  for 
multiplexing,  A/D,  timing,  and  control.  The  multiplexer  is  capable  of  data 
rates  in  excess  of  40,000  readings  per  second.  The  azimuth  error  due  to 
multiplexing  250  readings  per  revolution  on  three  channels  is  0.00675  degree, 
which  is  considered  nil. 

SHAFT  AZIMUTH  ENCODING 

The  data  recording  for  this  program  needs  to  be  referred  back  to  a  physical 
location  on  the  drive  shaft.  This  requirement  is  satisfied  by  the  fabrication 
of  a  special  Hughes  Helicopters  developed  electronic  circuit  that  encodes  the 
shaft  into  250  equally  spaced  segments  around  its  azimuth.  The  use  of 
250  units  is  an  arbitrary  but  convenient  buffer  size  in  the  digital  computer. 

The  encoding  system  consists  of  an  aluminum  disc  approximately  8  inches 
in  diameter  and  0.  065  inch  thick.  The  disc  is  machined  to  provide  250 
equally  spaced  radial  cuts  0.020  inch  wide  and  0.  25  inch  deep  from  the 
outside  diameter.  One  of  the  cuts  (the  master)  is  0.  50  inch  deep.  A  pair 
of  infrared  optical  encoders  arc  used  to  start  and  stop  a  pulse  circuit.  As 
the  master  slot  in  the  timing  plate  passes  the  optical  encoder,  an  automatic 
circuit  is  enabled,  providing  the  computer  with  a  trigger  at  each  of  250 
equally  spaced  azimuth  locations.  The  automatic  trigger  circuit  is  dis¬ 
abled  after  one  complete  revolution,  and  no  additional  readings  can  be  taken 
until  the  reset  switch  is  activated.  Figure  15  shows  the  encoding  disc  and 
th  e  operator  control  switch. 
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At  each  trigger,  the  computer  scans  the  required  data  channels,  stores  the 
information,  and  waits  for  the  next  trigger.  Appendix  A  provides  a  complete 
explanation  and  schematic  of  this  system. 

DIGITAL  COMPUTER 

The  HP  1000  digital  computer  is  programmed  to  do  the  following  tasks: 

a.  Calibrate  the  incoming  voltages  and  convert  the  readings  into 
engineering  units. 

b.  Read  the  data  during  a  measurement  run,  place  it  into  an  array  in 
memory,  and  print  the  results  on  a  line  printer  after  each  azimuth 
scan  (1  per  station  for  100  stations). 

c.  Graph  the  results  of  a  complete  run  on  a  plotter  and  print  the 
required  balance  correct,  ons  in  tabular  form  on  a  line  printer. 

Appendix  B  provides  listings  and  explanations  of  all  programs  used  for  this 
project.  Programs  are  written  in  both  BASIC  and  FORTRAN.  BASIC  is  used 
when  speed  is  not  a  factor,  while  FORTRAN  is  needed  to  take  data  and  per¬ 
form  the  computations  quickly. 
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SYSTEM  CALIBRATION/ RESOLUTION 


The  initial  run  of  this  system  was  checked  against  known  inputs  for  reliability, 
repeatability,  and  accuracy.  The  first  check  assured  that  the  voltages  read 
at  the  analog  output  of  the  Nova  Scope  were  being  recorded  by  the  computer. 
This  check  verified  that  magnitude,  polarity,  and  channel  designation  were 
correct. 

The  second  check  verified  that  all  the  program  logic  and  math  algorithms 
functioned  as  expected. 

The  third  check  measured  the  wall  thickness  and  eccentricity  of  a  calibration 
tube  (Figure  11)  made  specifically  for  this  program.  This  tube  was  measured 
with  micrometers  to  determine,  by  a  separate  method,  the  same  parameters 
that  the  ultrasonic  system  does. 

A  fourth  check  was  to  run  the  system  on  the  first  drive  shaft  specimen, 
S/N-1U,  and  keep  the  transducer  at  the  same  position  for  all  100  station 
points. 


This  test  proved  system  resolution  to  be  less  than  0.  001  inch  and  2  degrees 
for  locating  the  centroid.  The  results  of  this  test  are  plotted  in  Figure  4. 

In  an  ideal  system,  this  check  would  show  all  data  to  be  exactly  the  same. 


PROCEDURE  AND  RESULTS 


Two  Phase  1  AAH  shafts  (4.  5"  O.  D.  by  17  feet  long)  were  balanced  to 
demonstrate  the  technique.  The  first,  S/N-1U,  was  used  during  Phase  1 
development  of  the  YAH-64  in  the  ground  test  vehicle  and  was  not  flightworthy 
because  of  some  slight  damage  near  one  end,  caused  by  rubbing  on  a  frame. 
The  damage  is  slight,  a  photograph  would  not  show  it,  and  the  effect  on  the 
balance  process  is  considered  nil. 

All  balance  material  on  the  shaft  was  removed.  The  shaft  was  cleaned  of 
all  dirt  or  flakes  of  paint  or  epoxy  that  might  cause  an  ultrasonic  artifact 
and  the  shaft  was  installed  in  the  spin  fixture  for  the  ultrasonic  data  scan. 

The  first  scan  took  15  minutes,  and  produced  a  plot  and  tabulated  correction 
weight  chart  within  two  minutes.  The  scan  time  is  documented  by  the  data 
acquisition  rate  of  8  seconds  per  station. 

A  laboratory  technician  was  instructed  in  how  to  read  the  plotted  and  tabu¬ 
lated  data.  This  instruction  required  approximately  30  minutes. 

The  hand  lay-up  of  the  aluminum  balance  tape  used  as  the  correction  medium 
required  3.5  hours.  The  correction  process  is  accomplished  in  two  phases: 

a.  Lay  out,  using  a  suitable  marker,  the  azimuth  locations  (station  1 
to  station  100).  This  process  was  responsible  for  over  1.5  hours 
of  the  balancing  time. 

b.  Apply  the  tape  in  layers  and  lengths  at  the  required  azimuth,  start¬ 
ing  with  the  widest  tape  and  working  up  to  the  narrowest. 

The  initial  balance  was  not  good  enough,  and  upon  investigation  it  was  dis¬ 
covered  that  all  computer  corrections  were  based  on  0.  005-inch-thick 
aluminum  tape,  while  the  tape  being  used  was  0.  003  inch  thick.  A  subse¬ 
quent  run  of  the  same  data  with  a  new  constant  for  the  0.  003 -mch-thick- 
tape  provided  new  data.  A  new  taping  operation  was  done  in  a  1.  5 -hour  per 
period. 

The  operation  of  the  shaft  through  four  critical  speeds  (0-175  Hz)  was 
successful,  with  success  being  determined  by  dwelling  at  the  critical  speeds 
for  several  minutes  with  no  damper  smoking  or  evidence  of  excessive  heat¬ 
ing,  Figure  12  illustrates  the  shaft  and  mounting  features. 

The  results  of  S/N-5U  were  similar  to  S/N-1U  except  only  one  taping  was 
required. 
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The  adhesive  quality  of  the  aluminum  tape  used  in  this  demonstration  was 
great  enough  to  guarantee  adherence  to  5  times  the  maximum  operational 
rpm  of  the  YAH-64  tail  rotor  drive  shaft.  To  keep  the  adhesive  safe  from 
environmental  factors,  a  moisture-resistant  sealant  must  be  painted  over 
the  tape  and  especially  along  its  edges.  The  density  of  the  cured  epoxy 
sealant  is  0.  47  lb/cu.  in.  ,  and  has  been  found  to  have  little  or  no  effect  on 
the  balance  when  used  in  the  quantities  required  here. 

Table  1  presents  results  of  measured  lateral  response  at  the  aft  end  of  the 
shaft  bearing  mount.  Data  was  recorded  as  velocity  and  converted  to  dis¬ 
placement  and  acceleration  by  analysis.  This  table  is  presented  to  show, 
quantitatively,  the  dynamic  response  characteristics  of  the  drive  shaft  as 
mounted  in  the  spin  fixture.  Comparisons  to  response  characteristics 
when  used  in  another  environment,  such  as  an  A /  C,  may  show  significant 
differences  due  to  the  dynamics  of  the  new  structure  attached. 

Table  2  and  Figures  5  and  6  are  computer  generated  charts  and  graphs  of 
S/N  1U  and  5U  showing  mass  distribution  measurements  recorded,  and  the 
required  corrections.  Figures  7  and  8  are  spectrum  analysis  plots  of  before 
and  after  balance  corrections  on  S/N  5U.  Figures  9  and  10  are  spectrum 
analysis  plots  after  correction  runs  on  S/N  5U.  These  figures  show  the 
pertinent  vibratory  response  characteristics  measured  during  these  tests 
and  are  the  source  of  the  data  presented  in  Table  I.  The  effect  of  the 
balancing  process  is  clearly  shown  by  comparison  of  the  magnitudes  of  tin- 
various  plots. 

TABLE  l.  VIBRATION  CHARACTERISTICS  -  S/N-5U 


F  requency 
(Hz) 

Voltage* 

(rms) 

Velocity 

(in/sec) 

Displacement 

(Inches 

Peak -to- Peak) 

Accele  ration 
(g's) 

13 

0.  042 

0.  074 

0. 0018 

0.  016 

50 

0.  142 

0.  741 

0. 0047 

0.  603 

81** 

0. 0035 

0.  006 

0. 000024 

0.  008 

106 

0.  110 

0.  194 

0. 0006 

0.  3  34 

175 

0.  110 

0.  194 

0. 0004 

0.  550 

*Voltage  to  engineering  units  -  0.  00176 
**Ope  rational  speed  of  driveshaft  in  YAh 

in/sec  rms/mv  rms 

-  o4 
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l.  DIGITAL  BALANCE  PRINTOUT  -S/N5U 


V  *.  1980  II lit:  V  t:  1  11  S/N  5U 


1 


■4 


I  U 

l  I 
l 2 
1 

14 

r.. 


i  ■ 
i . 


i' 

•ij 


jit 


1  *yn  I  1111  :  I  <  PH 

l  Un  Ii  l  2  III*-  Ij  1  4  lipli  I  •  *  i  • 


l  ij  H 

1  it  u 

1  1 1  1/1 

1  41  i'i 

4  ii  ii 

ii  I  1 

t  it  ii 

I  ij  fi 

ij  l  ii 

ij  1  ii 

ii  I  I 

0  1  1 

0  I  1 

lj  II  1 

Ij  II  J 

lj  I  II 

lj  1  II 

lj  I  1 

‘I  I  I 

lj  I  II 

II  I  | 

ti  ;  i 

i  .i  .1 

i  ii  ii 

||  i  i 

it  i  i 

‘j  i  ii 

n  i  i 

i  ii 

i  "  i 

.  ii  i 

i  n  i 

i  ii  i 

i  i  i 

‘I  i 


i  ti  i  j  1 *  '  tie  a  i 


192 

189 

193 
191 
171 

1 9 1 
212 

190 
1  70 

144 

1 33 

148 

149 
121 

82 


*ji. 


i'i 


,4 

C-0 


36  1  ll 

37  i  : 

:3  l  I 

5 '5  >  II 

40  i 


1 

! 

I 

i.i 

ii 


89 


4 1  2  ‘I 

•f  1  1 

4  3  2 

4  4  2  'I 

49  2  U 


80 


81 


i.i 

ii 


ECONOMIC  ANAL\  S ilS 


A  cost  effectiveness  analysis  was  conducted  to  determine  the  impact  of  using 
this  method  of  balancing  the  YAH-64  supercritical  tail  rotor  drive  shafts. 

The  cost  to  install  the  system,  including  computer  and  programming  time, 
is  estimated  at  $51,  500.  The  estimate  is  in  1980  dollars  and  is  based  on  the 


following  equipment. 

Manufacture  r 

Model 

Name 

Cost 

HP 

9  82  5  A 

Desk  Top  Computer 

$12,  000 

HP 

3437A 

System  Voltmeter 

2,  000 

HP 

Typical 

Printer  Plotter 

5,  000 

HP 

HP -IB 

Interface  Card 

500 

HP 

A/R 

I/O  Roms 

1, 000 

HP 

3439  A 

Scanne r 

2,  000 

NUT 

3  000 

Nova  Scope 

18, 000 

Software  Development 

10, 000 

$51, 500 
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The  cost  of  this  equipment,  amortized  over  1,000  parts,  is  $51.50  each. 

The  cost  could  be  reduced  considerably  if  the  operation  were  installed  at  a 
facility  that  could  time-share  a  computer.  The  actual  C.  P.  U.  time  required 
to  do  this  job  is  very  short;  however,  during  the  data  scan  at  a  given  station, 
the  computer  must  devote  continuous  time  to  data  acquisition. 

The  minimum  equipment  required  is  the  ultrasonic  system,  the  software,  a 
high  speed  multiplexer,  and  the  A/D  converter.  It  is  estimated  that  this 
minimum  system  could  be  installed  for  approximately  $35,000. 

The  time  to  actually  balance  a  shaft  can  be  divided  into  two  areas:  data 
acquisition  and  balance  correction.  The  data  acquisition  and  printout  has 
been  demonstrated  to  take  less  than  15  minutes.  The  balance  correction 
has  been  demonstrated  to  take  approximately  3  hours.  These  are  actual 
documented  times  from  the  work  done  on  this  program;  however,  improve¬ 
ments  of  at  least  2:1  can  be  realized  under  production  conditions. 

For  comparison,  the  time  to  balance  the  YAH-64  tail  rotor  drive  shaft  using 
current  technology  is  documented  in  program  records  at  8  hours  per  shaft. 

Based  on  1000  shafts,  the  labor  saved  is  at  least  3.  75  hours  per  shaft.  As¬ 
suming  the  value  of  labor  through  G  &  A  to  be  $30.  00  per  hour,  this  repre¬ 
sents  a  total  labor  savings  of  $112.  50  per  shaft,  or  $112,  500.  00  for  1000  shaft 

The  YAH-64  uses  two  supercritical  and  two  subcritical  shafts  to  which  this 
process  could  be  applied.  Since  four  shafts  per  A/ C  represent  four  times 
as  much  balancing  labor,  as  mentioned  above,  the  costs  for  balancing 
1000  shipsets  of  Tail  Rotor  Drive  Shafts  could  be  reduced  approximately 
$450,  000.  00,  using  this  method. 
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CONCLUSIONS 


It  can  be  concluded  after  reviewing  the  data  generated  during  this  program 
that: 

a.  The  method  described  and  demonstrated  herein  does  locate  the 
center  of  gravity  of  the  shaft  on  the  center  of  rotation,  and  it  con¬ 
stitutes  a  condition  of  balance  which  is  not  changed  by  rotational 
speed. 

b.  The  degree  of  balance  required  is  dependent  on  the  desired  opera- 
ti->nal  speed  of  the  end  use  of  the  shaft. 

c.  A  shaft  designed  to  operate  supercritically  can  be  balanced  by 
making  one  quasi-static  data  run.  There  is  no  requirement  for 
high  speed  spin  runs. 
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RECOMMENDATIONS 


It  is  recommended  tnat  a  study  be  conducted  to  determine  the  practical 
limitations  of  shaft  interchangeability.  Although  a  shaft  can  be  balanced  to 
locate  its  center  of  gravity  within  0.  001  inch  of  the  center  rotation,  in  a 
balance  machine,  the  balance  is  only  as  good  as  the  location  of  the  center 
of  gravity  in  its  final  use.  This  should  be  done  during  the  original  design 
of  the  shaft  application,  so  that  balance  operations  can  be  done  on  a  bench 
set  up  during  manufacture.  The  eccentricity  of  all  mating  rotating  parts  in 
the  end  application  can  sum  to  a  significant  unbalance  which  can  only  be 
accounted  for  in  an  in-place  balance. 

This  program  has  highlighted  the  value  of  this  technique  as  a  receiving 
inspection  tool,  as  well  as  a  balancing  technique.  It  is  also  recommended 
that  this  system  be  considered  as  a  method  of  determining  the  acceptability 
of  any  tubular  materials  which  are  used  in  applications  where  straightness 
and  wall  thickness  uniformity  are  essential. 

Further  development  is  recommended  to  extend  the  automation  process, 
first  to  a  computer-driven  printer  that  locates  the  centerline  of  the  balance 
tape  along  the  drive  shaft,  and  in  a  second  step  to  develop  an  automatic 
tape-laying  mechanism. 
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PREAMPUF 


Lgure  1.  Schematic  of  balancing  system. 


Figure  2.  Schematic  of  balancing  process. 


Figure  3.  Ultrasonic  signal. 


STATION  NUMBER 


Figure  6.  Computer  plot  of  shaft  mass  distribution  S/N  5U 
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appendix  a 

TAIL  ROTOR  DRIVE  SHAFT  POSITION  ENCODER 


Circuit  Operation 

The  following  description  of  the  position  encoder  is  referenced  to  the  wiring 
schematic  in  Figure  A-l. 

Manual  Mode 


Upon  power-up,  a  manual  reset  must  be  issued  to  set  all  registers  to  a  ready 
state.  The  first  1/rev  timing  pulse  (high)  produces  a  corresponding  high 
pulse  at  pin  8  of  74LS13(B);  simultaneously,  position  0  of  the  250/rev  produces 
a  high  at  pin  6  of  74LS13(A).  These  simultaneous  pulses  cause  register 
74279(A)  to  become  set,  providing  a  "high"  to  pin  8  of  7408(C),  and  register 
74279(B)  to  be  reset  which  locks  out  all  further  set  signals  to  register 
74279(A)  until  a  system  reset  is  initiated.  Each  successive  250/rev  pulse  is 
then  "added"  with  the  set  of  register  74279(A)  and  fed  to  the  transmission  line 
where  it  is  used  to  start  and  stop  the  HP  computer  pacer  (H-run,  L.-stop). 

The  next  simultaneous  250/rev  and  1/rev  pulse  will  then  clock  the  7474(A) 
into  a  reset  state  which  in  turn  resets  register  74279(A),  thus  inhibiting  all 
further  pacer  pulses  until  a  system  reset  is  initiated. 

Automatic  Mode 


Automatic  operation  is  basically  the  same  as  manual  operation  except  a  hand¬ 
shake  operation  between  the  encoding  circuit  and  the  computer  is  instituted  and 
the  "data  in  error  automatic  shutdown"  circuit  is  enabled.  In  this  case,  the 
pacer  S/S  signal  is  disconnected  from  the  transmission  line  and  is  first  used 
to  trigger  the  74121  one  shot,  which  in  turn  resets  counter  7490  and  register 
7474(B).  Secondly,  the  pacer  S/S  signal  is  "added"  with  the  negation  of 
register  7474(B)  to  form  a  new  pacer  S/S  signal  which  is  fed  to  the  trans¬ 
mission  line.  When  the  pacer  is  started,  pacer  pulses  are  received  via  the 
transmission  where  they  are  counted  by  the  7490  counter. 

After  five  pace  pulses  are  received,  the  output  of  the  counter  w  ill  go  high 
which  sets  register  7474(B),  inhibiting  the  pa<  er  S/ S  signal  and  stopping 
pacer.  The  pacer  then  remains  off  until  the  next  250/ rev  pulse  and  the 
sequence  is  repeated.  Initial  start-up  and  final  shutdown  are  the  same  in  the 
automatic  mode  as  in  the  manual  mode. 
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gure  A-l.  Drive  shaft  position  encoder. 


In  the  automatic  mode,  the  data  error  shutdown  circuit  is  also  enabled. 

Should  the  250/rev  timing  pulse  fall  low  before  five  pace  pulses  are  received, 
the  shaft  travel  will  be  considered  too  excessive  for  the  data  to  be  valid  and 
the  7402  logic  circuits  will  place  the  encoder  circuit  in  a  shutdown  state, 
thus  terminating  operation  until  a  system  reset  is  initiated.  Two  internal 
dip  switches  are  provided  to  switch  between  the  auto/manual  mode. 

Data  Rate 


The  data  rate  is  limited  by  the  response  characteristics  of  the  timing  sensors 
and  the  external  measurements  to  be  taken.  The  rest  of  the  circuit  is  capable 
of  operating  in  excess  of  1  MHz. 
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APPENDIX  B 


SUPER  CRITICAL  DRIVE  SHAFT  BALANCING  SOFTWARE 
_ SOFTWARE  DOCUMENTATION  AND  LISTING 


The  Drive  Shaft  Software  package  consists  of  four  interrelated  programs. 

The  first  is  the  Calibration  program  calle4  DCAL.  This  program  allows 
recalculation  of  the  calibration  constants  which  are  required  for  the  proper 
running  of  the  programs  which  follow  it.  To  accomplish  this  the  program 
requires  two  aluminum  plates  of  known  thickness  and  "Eccentricity" 

(actually  what  is  needed  is  the  difference  in  heights  of  the  two  plates  relative 
to  the  fluid  nozzle). 

After  recalculation  of  the  calibration  constants,  the  program  called  DSHAFT 
is  run.  This  program  reads  the  calibration  values  from  the  disk  and  lists 
them  on  the  CRT  terminal.  It  then  calls  the  two  programs  which  actually 
accomplish  the  task  of  determining  the  amount  of  aluminum  tape  required  to 
balance  the  shaft.  These  two  programs  are  called  DS1  and  DS2.  The  first, 
DS1,  reads  in  data  from  the  Nova  Scopes  and  calculates  the  cross-sectional 
area  required  to  balance  each  of  the  100  stations  along  the  length  of  the  shaft. 
The  program  then  returns  control  to  DSHAFT,  which  calls  DS2.  DS2  then 
determines  the  number  of  layers  of  each  width  tape  needed  to  balance  each 
station.  Once  these  values  are  determined  the  program  both  graphs  and 
prints  the  calculated  tape  information.  When  this  program  completes  it 
returns  control  to  DSHAFT  which  executes  an  END  statement  and  turns  over 
control  to  basic. 

DCAL 


After  calling  the  subroutine  NORM  to  normalize  the  2313  data  acquisition 
system,  the  program  prompts  the  operator  to  input  the  value  of  wall  thickness 
for  the  first  calibration  plate  and  the  relative  eccentricity  of  that  plate.  One 
of  the  two  plates  may  be  taken  to  have  zero  eccentricity  and  used  as  a  stan¬ 
dard  value.  The  second  plate  then  has  an  eccentricity  relative  to  the  first. 

The  program  then  waits  until  the  operator  is  sure  that  the  proper  calibration 
plate  is  in  place  under  the  fluid  nozzle.  When  the  operator  types  "G"  then  the 
program  begins  taking  two  sets  of  100  voltage  readings.  One  set  contains  the 
voltage  output  from  the  Nova  Scope  for  Eccentricity  while  the  other  contains 
the  values  for  Wall  Thickness.  The  100  readings  of  each  set  are  then  averaged 
to  produce  one  voltage  reading  for  Wall  Thickness  and  one  voltage  reading  for 
Eccentricity.  The  program  then  notifies  the  operator  of  any  voltage  errors 


detected  and  of  the  completion  of  the  first  reading  sequence.  The  operator  is 
then  asked  to  enter  a  new  set  of  calibration  values.  These  consist  of  the  Wall 
Thickness  and  Eccentricity  for  the  second  plate.  The  program  then  waits  until 
the  operator  types  "G"  before  taking  the  next  set  of  readings.  This  allows 
time  to  replace  the  first  calibration  plate  by  the  second.  Three  sets  of  100 
readings  are  then  taken  and  averaged  to  produce  three  voltage  readings. 

These  are  the  two  Nova  Scope  readings  plus  a  temperature  reading.  This 
third  reading,  the  temperature,  is  used  to  account  for  variations  in  the  read¬ 
ings  due  to  the  change  in  speed  with  temperature  of  the  reflected  signal. 

Using  the  two  sets  of  readings  and  calibration  inputs,  the  program  then  cal¬ 
culates  the  slope  and  Zero  voltage  offset  for  both  Wall  Thickness  and 
Eccentricity.  These  slopes  have  units  of  Volts  per  Inch  (Volts/Inch)  and  are 
calculated  using  the  equation  for  a  line:  Y  =  m*X  +  b,  where  "m"  is  the 
slope  and  "b"  is  the  offset.  These  two  parameters  are  computed  using  the 
following  two  equations  which  also  include  temperature  offset: 

Slope:  m  =  (input  plate2  -  input  plate  1 )/ (volts  plate2  -  volts  plate  1 ) 

offset:  b  =  (input  platel)  -  (m*(volts  plate  1 )  )  -  (temp.  volts)/100 

These  values  are  then  printed  on  the  external  printer  and  the  values  saved  on 
the  disk  in  file  DSDATA.  The  voltage  and  input  data  are  also  stored  along 
6  with  the  calculated  slopes  and  offsets. 

F 

Note:  The  sign  of  the  slopes  is  unimportant,  as  the  program  DSHAFT  sets  the 
sign  of  the  slope  to  correspond  with  the  observed  behavior  of  the  output  of  the 
Nova  Scopes. 

DSHAFT 


This  program  begins  by  reading  the  calibration  values  calculated  by  the 
program  DCAL  from  the  disk  file  DSDATA  and  displaying  the  values  to  the 
screen.  However,  prior  to  displaying  the  data  the  program  corrects  for  the 
sign  of  the  slope  values  by  setting  the  Wall  Thickness  slope  to  a  positive  value 
and  the  Eccentricity  slope  to  a  negative  value.  The  program  then  waits  for 
the  operator  to  hit  Return  before  calling  the  FORTRAN  subroutine  DS1.  When 
DS1  is  called  it  is  passed  the  calibration  data  from  DSHAFT  and  returns  the 
calculated  values  for  Radius  of  centroid,  Area  required  to  balance  the  shaft, 
Angular  location  of  the  Centroid,  Radius  opposite  Centroid,  and  First  moment 
of  inertia  of  the  centroid.  These  calculated  values  are  then  stored  on  disk 
file  DSDATA  by  DSHAFT.  Next,  they  are  passed  to  the  FORTRAN  subroutine, 
DS2,  which  completes  the  calculations.  When  DS2  completes  and  returns 
control  to  DSHAFT  the  program  cycle  is  complete  and  DSHAFT  ends. 
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DS1 


This  is  the  key  program  in  the  sequence.  It  is  responsible  for  computing  the 
location  and  magnitude  of  the  imbalance  at  each  station.  To  do  this  the  pro¬ 
gram  requires  two  calibration  slopes:  the  wall  thickness  offset  value  at  zero 
volts  and  the  Radius  of  the  shaft  at  zero  degrees  azimuth.  Using  these 
values  the  program  takes  2  50  readings  at  each  station  at  even  angular  divi¬ 
sions  around  the  shaft.  The  readings  are  triggered  by  a  change  in  voltage  on 
a  control  channel  which  acts  as  a  pacer  for  the  readings.  When  the  voltage  on 
the  control  channel  goes  from  a  value  below  2  volts  to  a  value  in  excess  of 
6  volts  the  program  takes  one  reading  on  each  of  three  channels.  The  first 
channel  is  the  Eccentricity  voltage,  the  second  is  the  Wall  Thickne s s  voltage 
and  the  third  is  the  Temperature  voltage.  The  program  then  waits  for  the  next 
low  to  high  transition  before  taking  an  additional  reading.  This  continues  for 
250  readings  around  the  shaft.  (The  program  acts  as  a  positive  edge  trigger 
device  for  data  acquisition.)  If  after  250  readings  it  is  determined  that  addi¬ 
tional  triggers  are  occurring  on  the  pacer  channel  then  an  error  condition 
occurs  and  the  program  prepares  to  retake  the  last  set  of  readings  after  the 
operator  hits  the  return  key. 

After  the  250  sets  of  three  readings  have  been  taken  the  program  calculates 
the  location  of  the  centroid  and  the  centroid  magnitude  as  follows: 

6 
B 


Area  =  (Wall  thickness)*(Radius)*(Delta  Theta) 

This  is  a  rectangular  approximation  to  the  area.  The  radius  used  is 
the  radius  determined  by  the  er.  entricity  minus  half  the  wall 
thickness. 

b.  The  centroid  is  determined  in  Cartesian  coordinates  using  the  sine 
and  cosine  to  determine  the  Y  and  X  components  of  the  moment 
vectors.  These  values  are  summed  around  the  circle  to  determine 
the  Centroid  of  the  entire  circular  slice  of  the  shaft.  These  values 
are  then  converted  into  Polar  coordinates  using  R=SQRT(  (X*X)  +  (Y*Y)  ) 
and  Theta=  Arctan( Y/X). 


a.  First,  the  voltage  reading  for  Eccentricity  is  converted  into  a 

Radius  value  and  the  voltage  for  Wall  Thickness  is  converted  into 
inches.  Both  these  values  have  a  correction  for  temperature  varia¬ 
tions  by  adding  the  temperature  voltage  divided  by  100  to  the  inches 
value.  (Note  this  is  an  empirical  correction.  )  The  centroid  is  then 
calculated  by  summing  the  individual  moments  about  the  circle.  The 
individual  moments  are  determined  to  be  the  area  of  that  slice  times 
the  radius  to  the  center  of  that  area.  The  area  is  determined  as 


4  ’ 


c. 


These  values  are  then  used  to  determine  the  amount  of  area  required 
to  balance  the  shaft  at  that  station.  The  tape  must  be  added  at 
180  degrees  from  the  centroid  of  the  shaft.  This  angle  is  determined 
and  the  actual  radius  of  that  location  is  calculated.  This  value  is 
later  used  to  determine  the  moment  of  a  piece  of  tape  at  that  location. 
The  Radius  of  the  centroid  times  the  total  cross-sectional  area  deter¬ 
mines  the  magnitude  of  the  imbalance.  An  equivalent  moment  must 
be  added  in  the  opposite  direction  in  order  to  balance  the  shaft.  The 
area  of  tape  times  its  moment  arm  is  this  balance.  Therefore,  the 
area  of  tape  needed  =  (Radius  of  the  Centroid)*(Total  cross-sectional 
area)/Radius  opposite  the  Centroid). 


DS2 


This  FORTRAN  program  is  passed  the  information  calculated  by  DS1  and  com¬ 
putes  the  number  of  tape  layers  needed  to  balance  each  station.  It  then  graphs 
and  prints  this  data  along  with  a  graph  of  the  radius  of  the  centroid  and  the 
angle  opposite  the  centroid. 

The  determination  of  the  number  of  layers  of  each  tape  is  based  on  an  iterative 
process  for  the  first  size  tape  and  then  additional  layers  of  smaller  tape  as 
needed.  No  more  than  one  layer  of  the  smaller  layers  is  needed  at  any  station 
due  to  a  binary  approximation  method  because  each  smaller  layer  of  tape  is 
1/2  the  size  of  the  preceding  layer.  Therefore,  only  the  tape  widths  1  inch, 

1/2  inch  and  1/4  inch  are  needed. 

To  determine  the  effect  of  each  layer  of  tape,  the  radius  of  the  shaft  at  the 
balance  angle  is  used  and  updated  as  additional  layers  of  tape  are  added. 


PROGRAM  DC AL 


BCAL  T =00004  IS  ON  CPI 951 3  USING  00012  BUS  R-001B 

0001  10  (•'EH  THIS  IS  THE  DRIVE  SHflF  T  BALANCING  CRL 1  BRilT  I  OH  PPOGRuM  "  DC  HL 

0002  20F ILE SO >  DSDflT ft 

0003  30 DIM  IC  2  3>WC  102  3,  EC  102  3,  TC  100  3 

0004  40RRINT 

0005  50C RLLPhTER < 1 » 

000'S  CARPI  NT 
0 0 02  7  0  CRLLN  0  R  H  ■ .  1 '» 

0000  80PPINT 

0000  90RRINT  “MHAT  IS  THE  MOLL  THICKNESS"; 

0010  100 INPUT  wen 

0011  1  1  OF'P  I  NT 

0012  120PPINT  ” MHftT  IS  THE  ECCENTRICITY"! 

0013  130 INPUT  ECU 

0014  140PRINT 

0015  1  SOP EM  THIS  GOSUB  WRITS  UNTIL  YOU  ARE  READY  fo  TRIE  DATA 

001b  1  G  0  G  0  S  U  B  8 0 0 

0017  l?uLET  ZC 1  1=0 

0010  1S0LET  2123*0 

00 10  l'-’OLET  ES=0 

0o20  200CRLLAISQV<  100- -343, EC  3  3« E 1 > 

0021  2 1 OCALLAISOV* 100, -344, HC 33-E2  • 

0022  2 2 OPEN  THIS  GOSUB  CHECKS  THE  ERROR  RETURN  VhR1iiBI.ES  Eh 

0  0  2  3  2  3  0  G  0  S  U  B  S  5  0 

■•■024  2401 F  E 1  +E2=0THEN  270 

0025  25  OF'P  I  NT  “CHECK  INPUTS  HND  RETAKE  FIRST  SET" 

R02V.  2G0G0T0  ISO 

0027  270FRINT 

0027:  2.-30  PR  I  NT  ”  FIRST  READING  SET  DONE" 

0  7: 2"*  2  '  3  FOR  I  =  1 T  0  1 0  0 

On 20  SOOLET  ZC  1  3=2C  1  3  +  EC2+I  3 

00.:.  310LET  2C  2  3=ZC  2  3+Ht  2+ 13 

00  22  S.-.'ONtXT  I 

00  3  3  330LET  2C  1  3*ZC  1  I-  luO 

0024  J40LET  ZC  2  3*ZI  2  3  100 

0025  SSuPPINT  "NEXT  SET  OF  READINGS" 

0030  3S0RR1NT  "MHftT  IS  THE  WALL  THICKNESS" 5 
Ou37  270 INPUT  WC 21 

O'.1 37'  s.i'iER  I  NT  "WHRT  IS-  1  HE  ECCENTRICITY"; 

OoS4*  320  INPUT  CC  2  3 

..U4t.t  4C0RFM  THIS  GOSUB  WRITS  UNTIL  FuU  ARE  RLlll'Y  10  I  Hi  E  DATA 
>.i  ,i  4 1  4 1  ij  G  0  S  U 13  S  0  0 

0042  420CRLLH  '.SOVv  100,  -  24  El  3  1- E  1  1 
0043  4  3»i:RLLA;S0V-  100, -244, W!  ::1>S2.’ 

0044  i4uCRLLftiS0V-  100,  -045,  Tl  1 

0045  450PEM  THIS  GOSUB  CHE  Cl  3  THE  ERR  UR  RE  I  URN  VhR  i  RBI.  ES 

0u4b  40UG0SI3B  350 

0047  470 IF  E1+E2+E3=0THEN  500 

0040  4  7:  OF'P  I  NT  ”  CHECK  INPUTS  HMD  RETAKE  SECOND  SLC" 

00 42  4'RuGOTO  410 

0050  500PPINT  "DONE,  SRTOHTi  SFT  COMPIFTF." 

0051  510F0R  I =  2 T 0  100 

0052  520LET  EC  3 3=EC  3  I+EC  1  >2  I 
0053  SOOLET  WC  3  I=MC  3  3  +  WC  1+2  ] 

0054  540LET  Tl.  1J=TC  1  3+TI  I  3 
0055  550NEXT  I 
005b  5bOLET  EC3  3=EC3]  100 
0057  570LET  WC  0  3*WC  3  3-  100 
005S  5S0LET  TC I 3=TC 1 3  10000 

GU5?  5 2 OR EH  T  IS  THE  TEMP  COHSTftNT  lit  INCHES  CHANGE 
00 CO  EOOLET  M=-  WC  2  3 -HC  1  3  •  -’Ml  3  3-ZC  2  3  > 
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PROGRAM  DCAL  (CONT) 


00  tl 

00t-2 
00  £>.< 
00»54 
00t‘5 
0  0  *5 t" 
0  0 7 
006S 


o  :cu 

0071 


50:74 


i_.  7  7. 


007'!’'' 

C08u 

0031 

u.V;,,; 

7.1m;' 4 

7  USS 

000ft 

0.3:57 

00:53 

00:50 

0)300 

0091 


10LET  E  =  ■  EC  £  3 -EC  1  ]  1  CEL  3  Wl  1  J; 
20LET  E0  =  El 1 3-ZL 1 3-E-TC  1  ] 

30LET  WU=WL 1  3-EL  2  3>W-TC 1  ] 
40PPINT  #S 

THE  FOLLOW 1 Hi 


6P0PPIMT  *€■; 
GEOC  ALLI'ATER: 
870PRINT  ttE 
S80PRINT  #65 
690PRINT  #65 
700PR I  NT  St.! 
7 1  OFF  1  NT  HE: 
720PPINT  #E.; 
730PRIHT 
740  INPUT 


CONSTANTS  WE PE  ITU  LN  ON: 


"OK 

At 


WALL  TH I  Cl  NESS  CONSTANT 
WALL  TH ICINESS  OFFSET 
■ECCENTRICITY  CONSTANT 
'ECCENTRICITY  OFFSET 
'TEMPERATURE  OFFSET 
<  Y  o  r  N  >  « 


1  l  nc  1 1  •  :  V*  ■  1  *  1  -  "  5  W 

■  I nche  >  -  TWO 

*  I  f  t 'I  Tl *z  :  Vo  1 '  '  =  "  :  E 

■  Inc he i .■  -  “EO 

C  T €  tip-  1  '300  1  =  -;tc  n 


A! ="N" THEN  90 
NT  IS2<  1 


750  IF 
70 OPR  I 

77CPRINT  #2:  E0>  WO*  E<  W<  WC  1  ].WE21>EL  1  J.EC21>TL  1  J»ZL  1  l»2C21.£t  5  MIC  53 
7S0CHAIN  "I'SHAFT" 

7 90 REM  WAIT  COSUB 


SOORPINT  "When  you  ore  ready  to  take  data  t.re  G  and  Ret  run. " 

310 INPUT  AT 

'320 IF  HI  #"G"THEN  310 

3 SORE TURK 

34uRr.l1  VOLTAGE  ERROR  SUBROUTINE 

856  IF  ElttOPRINT  “Volta?*  error  on  Eccentric  it.  «•  re  ad  in?  Channel  343 

860  IF  E2U0PRIHT  "Volta?*  error  on  Wall  thick  nest  read  in?  Ch.  #344" 

070 IF  ESSOPRINT  "Voltave  error  on  Tem»erat  ure  readme  Channel  545 

S30PRINT  "VOLTAGE  ERROR  CHECKSUM  =  “5EH-E2+E3 

S90PR I  NT 

SCORE TURN 

910EHD 
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PR OGR  AM  DSHAFT 


DSHAFT  T =00004  IS  ON  CR19513  USING  06012  BLKS  P=0nlfc 

6001  16REM  THIS  IS  THE  DRIVE  SHAFT  BALANCING  PROGRAM  "DSHAFT" 

0002  20F I LESO > BSBAT  A 

0003  30REM  lit  *  1=BHGHTC  *].AC*  3=ATH£T At  *  3  -  PC  ^  3=P'£;Ak ■  E  ECO »  B=BAL 
0004  40D1M  Et  100  3-  6t  106  3- Wt  loo  3-  Rt  100  3-  At  lOu  1 

0005  50READ  #2-  1 
0  0  0  6  6 0 FOP  I  =  1 T  0  1 3 

0007  70PEAD  #2-ECI3 

0003  sonext  I 
0009  90PEH 

0016  106PEM  THE  FOLLOWING  TWO  ASSIGNMENTS  hPE  rn  ACCOUNT  FOR  THE 
0011  1 1  OPEN  FACT  THAT  THE  SLOPE  OF  (HE  ECCENIPl*  ITV  VOLTS  SHOULD  OE 

0012  1 2 OREM  NEGATIVE  WHILE  THE  SLOPE  OF  THE  MALI  (HICKNESS  VOLTAGE 


00 1 3 

130REN  S 

HOULD  HAVE  H  POSITIVE  SLOPE 

0014 

140LET  EC 

3  3=-ABS*  EC  31.- 

0015 

150LET  EC  4  ]=ABS<EC  4  I  > 

0016 

1 6Ur  P I  NT 

"OF SET  VALUE  FOR  ECCENTRIC  IT  Y-=  "Ifc'lll 

0017 

170PRINT 

"ECCENTRICITY  SLOPE*  In  V  •«  “  *  Et  3  J 

0013 

ISOPRINT 

"OF SET  VALUE  FOR  WALL  THICM(F.SS=  " :  1  t  81 

0019 

190PPINT 

"WALL  THICKNESS  SLOPE  *  In  V  •=  "!f.C41 

0020 

200PPINT 

"  CALIBRATION  VALUES  USED" 

0021 

21  OPR  I NT 

0022 

220PRIHT 

"WALL  THICKNESSES  USED  t  FIRST  =  "  !  Et  5  3 

"  SECOND 

0023 

230PRINT 

"ECCENTRICITIES  USED  s  FIRST  =  " : Et  7  3 

i"  SECOND 

0024 

240PPINT 

"VOLTAGES  READ  :  FIRST  SET  SECOND 

SET" 

0025 

250PRINT 

0326 

260 PP I  NT 

"CHANNEL  343"  (TAB*  13.*:  Et  10  3!  TAB*  :r.  •  5  lit  12 

3 

0027 

270PP I  NT 

"CHANNEL  344" !  TAB*.  13> !  EC  1 13!  1 AB*  35  •  5  Et  1  3 

3 

0023 

230PPINT 

0029 

2 90 PR  I  NT 

"TEMPERATURE  OFFSET" i EC  9  3 

0030 

300PP I  NT 

00  31 

31  OP PINT 

00  32 

320PRINT 

"WHEN  READY  HIT  RETURN" 

0033 

3.30  INPUT 

AT 

0034  340CALLBATER* 1  • 

0035  350C ALLBS1 *  EC  1  3-Bt  l  3  -  WC  1  3»  PC  1  ].  At  1  3* 

0036  3G0C  ALLDATEP 1 6  1 

0037  370PEAD  #2  >2 

0033  330F0P  1  =  1  TO  100 

0039  390 IF  I =51  READ  #2*3 

0040  400PPINT  #2-ECI3 

0041  4 1  ONE XT  I 

0042  420PEAB  #2-4 

0043  4  30F0R  1  =  1  TO  100 

0044  440 IF  I =51  READ  #2-5 

0045  450PR  i  nt  #2-em 

0046  460ne:;t  i 

0047  470PEAD  #2-6 
0043  4S0F0R  I  =  1 T 0  100 
0049  490 IF  I=51READ  #2-7 

0050  50OPRINT  #2!WCI3 
0051  51  ONE XT  I 

0052  520READ  #2-8 
0053  5 30 FOP  1  =  1  TO  100 

0054  540 IF  I =51  READ  #2-9 

0055  550PPINT  #2 5  PC  II 
>305*6  560NEXT  I 
0057  570PEAD  #2-10 
005:?  5S0F0P  1  =  1  TO  100 
0059  590 IF  I =51  READ  #2-11 

0060  60uPPINT  #2- AC  13 
0061  61  c-ne; ;t  i 

0062  620PA0St 

0063  G3uC ALLBS2*  Et  1  I-Bt  1  3-MC  1  I.  PC  1  I- At  1  I  ■ 
0064  640ENB 


EC  6] 
EC  8  3 
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PROGRAM  DS1 


BS1  T =C'0tt04  IS  OH  CR 19513  USING  00038  BLKs  F=Oj3-t 


0001 
0002 
o  i\  0  3 
0004 
0005 

000b 

0007 
0000 
0 0 0 0 
00  i  0 
00  U 
0012 
0010 
oo  1 4 
0015 


•lV'1  $ 
y  :■  l  \« 
.  2  o 
0021 
0.,i.":2' 


•j  j 


0  o  2  ■  *■ 


j  1.1  >»*1 


r.  *. » * 


0040 
»"  04  1 

»c 

0  0  4 

0044 

?04S 
0*%4 
0  »J  4  "* 
0..4.; 

0040 
00c.y 
of»5 1 
0l»52 
005  ’• 
0  0  5  4 
0055 
0o5*’ 
0057 
005 


FTN4.L 

SUBROUT I  ME  DS 1 >  E C  L . BHL • BWGHT . REAR *  HlHtTO  • 

DIMENSION  VOLTS'  3*  250  .' *  V  75. '.ECO'  10»'  • .  CiiL  •  100  •  •  I T I  ■:  5  •  >  MT I M  <  5  > 

DIMENSION  BWGHT •  1 00 > *  RBAR *.  1 OO  •  >  ATHET A ■  100  ' 

DATA  PI  3. 141 502054 
DATA  THOP I  b . 28  3 1 85 3 
DOTH  DTHETA  .0251327- 
DAT  ft  ISTAT  1 

C  THIS  PROGRAM  CALCULATES  THE  CENTROID  OF  H  PLANAR  SLICE  OF  A  DRIVE 
C  SHAFT  TAKEN  NORMAL  TO  THE  AMIS  OF  ROIAHON. 

C 

C  THE  DATA  IS  TAt EN  IN  POLAR  COORDINATES  AND  THEN  CONVERTED  INTO 
C  CARTESIAN  FORM  USING  THU  RELA1 IONS: 

C  ::=rC  A';.  /  t  I  ,>1  <  '1  • 

C  Y  =  rS IN'* he »  o  • 

f  VOLTS* 1* I >  WILL  HOLD  THE  VALUES  FOR  E >  C  ■  I S T A 1  -LIITR I C I  TV 

f  V0LTS*2vIa.  WILL  HOLD  THE  VALUES  F OF'  HULL  THICtHESS 

C  VOLTS'- 3*  I  1  WILL  HOLD  THE  VALUES  FOR  WATER  1 EHT ERATURE 

C  RADIUS  WILL  HOLD  THE  LOCATION  OF  THE  Me  l  OR  THE  SECTOR 

C  AREA  WILL  HOLD  'I  HU  C  ALCULATED  AREA  OF  EACH  SEi.TuR 

C  THE  FOLLOWING  ASSUMPTIONS  ARE  MAUL: 

C  1.  THE  SECTOR  IS  APPOM.  A  RECTANGLE 

C  2.  THE  CENTER  OF  MASS*  ID'*  FOR  I  hi  -ax  TOR  IS  HALF 

C  THE  WALL  THICKNESS  FROM  I  HE  aU i UR  EDGE. 

C  3.  THE  SECTOR  HAS  UNIFORM  DENSITY 

C 

C  XA  IS  THE  SUMATIOH  VARIABLE  FOR  X  HULIIhLIED  BY  AREA 

C  VA  IS  THE  SUMATIOH  VARIABLE  FOR  V  MU!  V 1  PI  I  ED  BY  AREA 

C  A  IS  THE  SUMATIOH  VARIABLE  KuR  IHF  i iF'f  A 

r 

C  THIS  BEGINS  THE  DUMMY  ARRAY  VALAI  Sl.'  TIOH  FOR  DEBUG 

C  THIS  BEGINS  THE  INITIALIZATION  SE'.TIuii  1  UR  DEFINING  i.ONSTANTS 

WRITE' 1*2 ' 

READ  Cl**'*  1  PRINT 
WRITE <1 *3' 

READ  *  1  *  * >  RC ON 

RCON-RC'ON  2 
WRITE'  1* 13  ' 

READ  '1  *  *)  PF 

VI =ECC ■  l  ' 

CONI  =F.CC  •  3  • 

V2-ECC • 2  ■ 

C0N2  =  ECC ■ 4  • 

WRITE’  1*  17 •CONt.CON2.RCON 
PMA>:=6. 0 
PMIN*2 . 0 

C  THIS  IS  THE  DATA  ml OU I  .- 1  f  1 1 -M  L"OK 


CALL  NORM 

•  i  • 

call  e 

::EC 

-ii* 

I  T  l  *  l 

YEAR  ■ 

DO  50 

I  - 

•  %  ;  r 

n 

CALL 

hi  : 

or  -  i 

F'-.-'nl  r ,  | 

VI  FT  • 

i  F  'P 

VOLT . GT 

.('Mill 

f,  i.i  1  0 

•'ALL 

A  IS 

OF  ■  1 

«  ::4*:  « 

K'-'F'l  1  *  I 

"I  F  F  • 

IF  >P 

VOLT.LT 

.  F  ha: 

•  ' ,  0  T  ii 

r_. 

CALL 

HIS 

OF  -  ♦ 

3 .  34  . 

..VOLT  S' 

i*i* 

CONTINUE 

CALL 

F-hF 

EF  ■  1 

•  .H. 

1  ■  I  ERR 

CALL 

H  I  ■;* 

. F  •  - 

7k . .  - 

.4i-..*V.  It'FF-  ' 
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PROGRAM  DS1  (CONT) 


0060 

C  ALL  FACER c 0 -  0 > 0  •  1 .  I EPF'  '• 

00*51 

VHIN=V'  1  • 

0062 

vmax=v*  i  :• 

006  j 

DO  27  1=2-  7 S' 

006  J 

IF  <  V <  I  > .  LT .  VNIH  ■  VI1IN=V'  I  ■ 

Oil 

if  <: vc  i  '.gt.  viih: :  ■  vi-ifi::=v*  i  • 

0060 

27 

com  I  HUE 

0067 

IF  CCVMAX-VHIN  •  .LT.  2*  GOTm  2:: 

0  0  6  y 

WRITE'  1. 24..  IS  TAT 

0  *’*■  6  1  > 

F  E  H  D  •.  1  >  * )  HALT 

0070 

GOTO  1 

*  9  7  l 

C 

COLL  EXEC'.  1  1  .  tl Till.  i'.'EAR  • 

•''072 

c 

::min=<mtiih  3 .-I  r  i  ■  3  •  '*.20.0 

00  79 

c 

sec=xm  1  n+mt  1  i-i •  2  •  - 1  t  1  c  2  >  <  1  ■  ht  1  m-  1 

0074 

c 

WRITE  '  6 *  20  *  :-F C 

>2075 

r  3 

ET  SU MOTION  VARIABLES  TO  ZERO 

00  7  £ 

Kfi=0 

0077 

YA=0 

H  =  0 

0  i*i  7  9 

0000 

THIS  BEGINS  THE  CALCULATION  SEi  T ION 

oo  o  i 

c 

0032 

V 1  =  VOLTS'.  1.1  -CONI 

0083 

DO  ICO  1-1.250 

00S4 

I F  < F'F .  EC! .  1  >  liR I TE ' G « 1  ft  ■  I .  VOL  1  ■  1  .  I 

1,1 0 S 

VOLTS'.:  1  .  I  '-FCON  + VOLTS'  1  .  I  —rONl  VI  * 

0000 

V0LTS<2»  1  '-VOLTS' 2>  I  '*i'OII2*V2  ►  VOLT' 

0087 

RADIUS* VOLT S.'  1  •  I  ■  -  •  VOLT  S'  2.  1  .■  2  • 

OOo.. 

AF:£A=VOLTS '.: 2 *  I  -PhUIUS-DTHi.  It! 

00*9 

IF  CPF.  EC!.  1  >  WRITE' 0.  19  >  VOL  TS'  1 .  1  . ' 

0090 

THETA*'  1-1  •-■D  THETA 

OU'9  1 

::a=  xa+' arra-radiu'.-  cos- tmcta  ■  • 

0  0  9  J 

YA=  YA  ■ AREA * RADIUS *S IN'  THE ih  •  • 

0093 

A*  A+APEA 

in '*4.; 

100 

CONTINUE 

0  0  9 

xbar*  ::a  a 

0096 

YEAR*  YA  H 

0097 

RBARC I  ST  AT  »*SORT'.  •  ::bar*xbar  »  +  '  YBARVi'I 

0090 

THETA=AT  AN2 <  YBAR  > XBAR > 

0099 

c 

0 1 00 

c 

BEGIN  PRINTOUT  SECTION 

0 1 0 1 

T  H  E  T  A  =  T  H  E  T  A  *  3  6  0  T  W  0  P  I 

0102 

IF  'THETh  .LT.  O'  THETA* 3G0+ THETA 

0 1 0  j 

WRITE  Cl  PRINT. 7 •  IS  TAT 

0104 

WRITE'  I  PR  I  NT- 9  ■ 

0 1  OS 

WR  I  TE '  I  PR  I  NT .  1 0 ■  1  PAR •  1 STAT  . 

0 1 06 

WRITE'  IPRIin. 11  •  THETA 

0 1 0  7 

WRITEC I PR I NT. 12 ■  A 

0 1  0 ... 

IF  'THETh  .  L  r .  lftii'  ANGLE  1  META  *  1  :':n. 

0109 

IF  C THETA  . GE.  ISO-  ANGLE  =  THE  1 M - 1 : :A. 

0110 

THETA=ANGLE 

0111 

ANGLE*' ANGLE-25 •  30 

0112 

ITHETA* I  NT • ANGLE ■ 

0113 

DT =ANGLE-FI.  OAT  l  THETh  • 

0114 

ITHETA=ITHETh+1 

Oils 

I  ANGLE* l THETA* 1 

o  1 1  >; 

IF  '  I ANGLE  .  G'i.  450'  1ANG11  .1 

0117 

DELTA* VOL  Tv  1  •  i  ANGLE  '-VOI  1  |  .  |  i  lu'  ill 

0 1  1 .3 

EC  C '  1ST  AT  ■  =VOL  T'_.  •  1 » I  THETA  •  ♦ •  14  1  1  n  •  M 

0  1  1  4 

WF I  TE '  I P F I H T  •  14  •  T ML  X  M « t i  L  •  (  T  H  1  ■ 

0  1 

ATHETA'  1 2 1  n  T  . NET A 

0 1  2 1 

8hL*  I T li T  •  ~h  ♦  I  !  hF •  I  TmI  • 

0 ;  2  ; 

BRIGHT  •  1ST  A I  '=1  ML'  1  31  AT  ■  li  r  •  J  ■  >  |  n  '  . 

1 00 . 0 


.  VOI.  IS'  2.  I 
VOL  IS  ■  3 «  1  >  • 

1  3 «  l  '*.01 


.  Vo l 
.01 


•'ML  I  S'  2.  I  '  •  RADIUS .  MRL 1 1 


Pi  IK 
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PROGRAM  DS1  (CONT) 


0123  WRITE*  I  PRINT .  21  '  BNGHT  *!  13  IT)  i  • 

0124  IF  LABS*  REAP*  I  :-_-TAT  *  *  .LT.  O.uOol*  HR  111.*  I  PR  I II I  >  15'* 

0125  ISTAT=ISTAT+1 

0 1 26  WRITE*'  I  PRINT  >  16  * 

0127  IF  (ISTAT  .l.E.  100,*  GOTO  1 

0123  2  FORMAT* "WHAT  IS  THE  PRINTER  *1  uR  6>?"* 

0120  3  FORMAT*  "DIAMETER  IN  INCHES'5"*  * 

01  30  4  FORMAT*.  "Eccent  ricti  y  Constant  *  Ini  i,f=  Volt  ■*••">  • 

0131  7  FORMAT*  *  "Noll  Thi'il.  lotut.  *  lii'.hti  V  oil  .*'?"*  .* 

0132  6  FORMAT*  > "Wort r  TtnPt  rot  urt  * IK  -*t  •; ; .  V-*-*l  i  *''".:• 

0 1  3  3  7  FORMAT  *  " ST AT  1  ON  Nl  IHBER  " -  I  3 > 

0134  '3  FORMAT*  "CENTROID  IN  POLAR  COORDINATES* 

0135  10  FORMAT <“  RADIUS  IN  INCHES  =  4  > 

0136  11  FORMAT*"  DEGREES  A2IMNTH  * 

0137  12  FORMAT < Si; >  "TOTAL  N-SECTIOI llil.  mREh  -  "  .  F  10. 4> 

0130  1:3  FORMAT*.  "DO  YOU  WANT  A  DATA  DUMP  •  1  =Yl  S  * 

013’'-*  14  FORMAT <  "ADD  WEIGHT  AT  "»F5.2* "  DEG  A7IMUTH.  RADIUS  -• 

vl40  1  j  FORMAT'  "**  BALANCED  *■  * "  * 

0141  1*3  FORMAT*" _ _ 

0142  17  FORMAT C ECC  CONST*" * F8. 4  "  WALL  Cut  1ST  =  "  «  Fs  .  *1  "RAD  I  (JO  - "  •  E 

0143  18  FORMAT*:  I  3  "ECC  V.  .  FS.  5*  "  WALL  V.  •ES.  G,  "  TEMP  V.  "-E 

0144  13  FORMAT*.  "CCG *" . FG. 5- "  WALL  =  " >  F3. 5«  "  RADIUS*"  *  FA.  5*  "IIRI  M* 

0145  21  FORMAT*:  ”•  AREA**  M'.ENTRO  ID  RADIUS*  *  f  CCEN I  K 1 II Y"  . 

0146  t"  OPPOSITE  CENTROID)*"* FT. 5) 

0147  22  FORMAT*;  "WK  AT  IS  THE  OFFSET  VALUE  FOR  ECCENTRIC  I  TY”"  *  * 

014:3  23  FORMAT*. “WHAT  IS  THE  OFFSET  /'Hl.iH'  FOP  WA  L  THICKNESS"". 

0143  24  FORMAT <  "ERROR-—  MORE  THAN  250  I'UI  GEO  -Ei  lEVL.’V 

0150  T “RESET  FOR  STATION  ">I3.“  THEN  HIT  RETURN.",* 

0151  RETURN 

0152  END 

0153  I 
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PROGRAM  DS2 


I»Sj 

T = ft 0 0 0 4  IS  ON  C P 1  3 j 

IS  USING  OOOSS  Eli  .  P---0.. 

000 1 

FTM4.L 

0  0 0  4 

SUE: POUT  I  HE  In  4 

-  ec c < eal < eiigh i  •  fl hi-  .nriii 

0005 

iii mens; ou  i  ■  j  ■ 

.  I. ,  -  X  -  -  I V-  1  1  ■  •  ECl.  •  1 1'U  -  •  1 

0004 

DIMENSION  L  IGHT  1  oO  • .  FEf.F  -  1  On  -  >  Hi  Ml  1  n • 

000  c> 

DIMENSION  MIT- 

1  i  U.  I  ' 

0000 

INTEGER  TCT-  ION. S  -.El  I C . FT  I r 

o  o  o  r 

data  ema::  o.o 

0005 

DATA  P mm::  0.0 

000^ 

data  istat  : 

00 1 0 

data  i -  i  4 .  ah 

•-•T  -  SHUT  .  LHIO.  ,:IIN  . .  :HNU« . 

0011 

DRTR  I Y  jO*  4hDH»  jHH*. «  4 1 1 5  «  JUNE  « ,. 

Oo  1 J 

DATA  T  .3.. IS 

. .  0C25 

0015 

DO  2  1=1. loO 

OC  1 4 

TCT  I  -  1  -  =  0 

0  0 1  *?. 

TCT <  I  >  2  1  =0 

l  *  4  5 

TCT-  I  -  3  -  =  0 

i  7 

YHT •  I  -  =0 

OolO 

I F  (REAP  -  I  - .  G 

1  .  PMA'  :  ■  FMri,  API  iHh  -  1  - 

/ .  ’  *  •;« 

IF  -EHGHT- I • . 

gt.eha:;  -  bntlaisugnt-  i  - 

O.j  jo 

CONTINUE 

00  j  l 

1  3  FORMAT  GI  ■ 

;;  .'i  J  j 

‘  J 

•n 

T1 

3 

— i 

-I  l-'-n  1  iii-'.l’i  1  S  in-1 

>.'  % 

f  "  An?U  •.  Ii«=»r 

-  “  - 

00  j4 

lo  FORMAT  -  12X>  1  S< 

r::»  12.  12-0:;.  is*  14;:.  1. 

Oo  j‘*. 

B  T I C  =  1 0 

0  •.*  j  5 

ESC ALE  =  300 

k; .  j 

IF  iBKA.'.  AE. 

. 01  -  GOTO  sno 

l'v.  t’  » 

ESC  ALE  -  1  SO 

IF  (BMA::  .le. 

.  02  .-  GO  1  0  SuO 

i;-.i  ;.n 

B SCALE1 GO 

0  0  j:  1 

IF  -:  BMh: :  .  Lt. 

.OS-  GOTO  SO C> 

0  0 

etic=is 

0  0  j  /• 

ESCALE-40 

oo  ;.4 

if  -emh::  .le. 

. ors  -  GOTO  soo 

0  0 *•' 

.  j  o  *,  ‘ ' 

BT 1C  =  10 

E-SCALE-  so 

IF  1 EMNI.  .A. 

. 1-  GOTO  SOu 

Om 

ET I C  =  i  S 

00  • 

BSCALE-20 

0  ^  4  n 

IF  (  EMIT.:  .  LE. 

.IS-  GOT 0  SOu 

J'!U4  . 

£■•  T I C  =  1  u 

004.: 

BSC  ALE  -  IS 

004  j- 

300  PTIC=10 

0u44 

PSCALE =150 

0O4f. 

IF  -  PMA.  AE, 

.  .  02  -  GOTO  So -T 

0045 

PSCALE  =C0 

004  7 

if  crna:-:  ae, 

,  .OS-  GOTu  SOO 

0  0  4 

R  T I C  “  1 5 

0045 

RoCALE-40 

00  50 

IF  CRM  AM  .LE, 

,  .ors-  GOTO  SOO 

0  0  j  1 

PTIC-10 

0  n  S  j 

PSCALE -  o 

00c- 

if  -pith::  ae 

.  .  1  -  GOTO  Soo 

OOf  4 

P  T  I C  =  1 S 

0  0  A  ‘  • 

PSCALE -SO 

I'M'’  7 

IF  -PIUT.:  .lE 
PT I C  =  1  o 
PSCALE--  IS 

.  .  A-  GOTO  ji'O 

IF  -PMA.:  .LE 

.  .  -  '  Gu  1 1-1  Si  hi 

1  Hl_1  1 


i'Dii  In-  l  -ju> 


H  1 1  j 


.'HEP 
.jtIHT  - 


Jll-  h>  .MPL< .HA 


1  4  ln-'li 


1 


1 


1 


7 

B 


I 


PROGRAM  DS2  (CONT) 


\  M  H 

■:  ALL  PLTi  I.I'  ;  ■ 

C  ALL  LLKFf 

J  ; 

L  hll  r  hi  f  '  1  L_. .  •  1  h  .  1 

M.'lr  .j 

CALL  ri.C'T  ■  0  .  .  ti  . .  ■ 

*-  '  ’  • 

*:  All  plot •  c . 

C  PlL  C  .  Lit-  ■  R-  . 

; .  «  0 .  .  - . 

I'O  s  i c  [  =  i .  c . :  c 

i  ;  -  .  j  •  ■  >  .  c-  T  R  ■ 

CL  i-VMe-  .T-VR.I.l 

• .  «  n  .  «  - » 

‘  V. 

t- 

CONTihOE 

CALL  :•  V M i- 1  .  1  >  5 .  j •  ,  1 

11 .  « u .  .  .1 

I'O  5E0  1  =  1.'. 

i\  .  I  ^ 

: ;  = .  ^.+  •  1.4  - 1  • 

^i; 

call  cvmB'  i.f. .  i • 

|C  NT  I  if.  L 

•  .  *n.  *  -  j  • 

0* 

C  All  I-YMI.  ■  1-1.  ..  .  f. 

:  »  1  I  .  «  U .  « 

call  :•  it-it- •  .. 

i'C  f  1  -  i  >  E  1 1 

i  «  1  ;  .  »  0  .  * 

.  :  '  •  • 

VI  -  3.  c.  -  •  I  *  EC  1 1  1 

L  AH  'EY  HE:  ■  1  -  .  f.  V  1  .  0  . 

•  1  ;  .  .  r«M.  , 

•  : 

0  .•  . 

■  OUT  I  NOE 

i'0  j  ;  -  i , :  i_. 

1  =  14.0-.  i .  a  *  i  • 

CPU  :  i'HL  ■ . :  i  •  .  j  •  .  1  >  1 

: .  «  U  .  •  -.1  • 

•  >* 

c  ..  r. 

'"  OUT  I  HUE 
•  7  I C  =  1 

: :  I  =  .  c.4 

CALI  VMC-:  •  0 .  B..J  <  .  f  •  .  1  •  lj .  .  0 .  .  -  1  • 


I'O  f  oo  i  -- 1  • :  i 

*’  .  'Vi 

CALL  no  ml  ,:f . .  l . ,  r  .  o. .  - 1 

V  T I C  =  1-1 0 

r .  , . 

:i  =  '  1*1.4 ...  or*.. 

4  '  ■  , 

V-ir. 

''OUT  INI'!: 

•V.-4 

'-  ALL  VNL  ON  -  .  1  E  •  I  0 1  1  ■  .  i' . 

LHLL  ■:  VHt'  .  J  .  0 . 0 .  .  1  L  .  1  V '  1  .NO*. 

0  0 '?  c 

DO  fOO  1=0<4 

Cuj 

V I  = . £25+ •  1*0.0  ' 

000:.. 

'.'TIC*  I *0. 0  EEC ALL 

OuOO 

'  HLL  lilJMtv  .4.  VI,  .  1 .  YT  I C >  *1.1.  .  + 

0100 

C  *. 

■j 

CONTINUE 

0 1 0 1 

CALL  PLOT'  .  f .  5.  L  . 

0 1 OE 

DO  f 40  I=0-PTIL 

0 1 0  ;. 

VI  =  ;-:.f*'  t-j.  PTIC  ■ 

0104 

CALL  EYHB '  .  j  •  Y  I  •  .  1 ,  1  .  Afi.  ,  - v  . 

0 1  Of. 

-•40 

CONT I  NUB 

0100 

CALL  YHC '  .  1  f .  .  1  ■  1  1 .  .  ;t.  .  -E  ■ 

0 1 0  7 

DO  f f A  I  -ii.  in 

0100 

::i  =  .f+'  i*i.4  ■ 

0 1 00 

CALL.  AVNE..  ::  I  •  t .  f  <  .  1 .  1  5.  .  0.  ,  -E'  - 

0 1  1 0 

c,c.u 

CONTINUE 

0 11 1 

LhLL  :>VMP-'  1 4 . 0 ,  t* .  f  *  .  1 ,  1  1 .  .  0.  ,  - 

0 1  R 

DO  f BO  I =0- PTIC 

0 1  1  0 

V I  =  B . f - ■  I  ♦  i .  FT  1C  • 

0114 

CALL  EYMB 1  1 4 .  f  •  V I  •  .  1  >  1  i .  .  -ft.  ,  -E 

0  1  1  N 

J  K>  0 

C  OUT  I  hot. 

0 1  1  c 

CALL  MU'lhi.  V.  VEf.  V  V  .  1  .fO.  .ft.  .  - 

o  1 1  r 

CALL  NUMB'  .  •  J.  Of.  .  1  >  0.  >  00.  .  -  1  ■ 

.,  1  1  ■: 

do  fro  i=i>4 

0 1  1 

T  I L  =  i  »  .  B.  F  LhLL 

0100 

V  I  =  i.  E'Ef  *"  I  *  .  >:  ' 

•J  t  w  1 

CALL  NONE  ■  .4.  VI  ..  1  .VI  p  .•-it.  .  . 

0  1  i . 

j  r  y J 

C  OUT  I  Ml  .f: 
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PROGRAM  DS2  (CONT) 


o  i  2 : 

!  V  •  1 =  1 3 

0124 

IY-.  2>=2HPR 

0125 

IV.  3  ■  =2ri['  I 

0 1  26 

I <  4  ■  =2HU3 

0127 

IY‘..5>=2H  0 

0120 

IV  *.6>  =  3HF 

0120 

I  i’.  7>  =2HCE 

0 1  BO 

IV' 3  >  =2HHT 

0 1  3  1 

I Y- 0)=2HR0 

0122 

IY<10>=2HID 

,*1  l  ' 

C  ALL  SVMB'. .  2»  3. 0  -  .  1 

2-  IT  '  1  ' -  00 . •  1  • 

0 1  34 

CALL  F'LC'T'  •  5j  6. 5»  3  • 

0125 

lift  530  1  =  1-3 

0136 

YI=6.5  +  '  I  +  . 375  ■ 

0 1 37 

CALL  SVMB'. .  5  -  VI  -  .  1  - 

13. -00. --2  - 

0133  530 

CONTINUE 

0130 

Bn  500  1=0-10 

0140 

::i=.  5+  ■  i*i. 4> 

014  1 

CALL  SVMB C<  1-0.5- . 1 

-  13. -0. - -2  ■ 

0142  500 

CONTINUE 

0143 

BO  60 Ci  1=0-3 

0144 

VI=0.  5- I*. 375  1 

0145 

CALL  SYMB< 14.5- VI  - . 

1 . 13. > OO. - -2> 

0146  600 

CONTINUE 

0147 

I V  (.  1  '■  =  1 4 

0143 

IY<2>*2HBR 

0140 

I  Y<3>=2'MLR 

0150 

IYU>»2HNC 

0151 

IV '1 5  ,=2HE 

0152 

IT'*:  6  ’  =  2HhN 

0153 

I V  ( 7  ;*=2HgL 

0154 

IT'-IS  :-  =  2HE 

0155 

CRLL  NUMB -.7.  375-6. 

3-  .  1  -  50.  -  0.  •  •  1  > 

0156 

CALL  NUMB- .4-6.55- , 

.i-o.i oo . « -  r- 

0157 

BO  610  1=1-4 

015.3 

YTIC=OuM 

0150 

V I  =6. 33+ -  I *.75' 

0 1 60 

CRLL  NUMB* . 4 -  V I  -  .  1 

•  V T I C >  00 . > -  1  ' 

0 161  610 

CONTINUE 

0 1 6  2 

CRLL  SVMB'  .2-7.  . 

12-  IV-  1  '-00.  -  1  • 

0 1 6  3 

ChLL  PLOT' .5- .5- 

-  3  • 

0164  C 

THE  REMAINING  CnLl. 

3  POP  PLOTTING  BATA  Hli  1 

0165 

BO  730  1=1-3 

0 1 6  3  700 

C  OUNT  =  0 

0  1 6  7" 

I  STRT=  1 

►v  i  v  L- 

CALL  PLOT'  .  14-  VHT.  1  ■ .  3  ■ 

at no 

EK=ECC'  ISTRT  '  +  .U015 

0 1 70 

AM=EK*tX-»S  IN'  i  I  • 

E  V  ■  100 

0 1 7  i 

RUM  I N=Ei A E  V-  3  IN'  "1 

•3  -  i ■  uu 

0 1 7  7 

IF  <-.AX-R;:mn  ■  .c 

T.  BRL •  ISTRT  ■  •  GO  TO  -'-‘i 

0  i 

IF  •:  I  .EC.  1-  CO-JHT  =  1 

0174 

TCT'  ISTRT- I  -  =  rcr- 

I  ‘  *T  hT-  I  1  +  1 

7,  *  —  t 

VHT' ISTRT '=VHT'  I 

•‘  T at  ■  <  ■  o::*r  •  hi. e  es  ■ 

v  i  7 1 

EAL I '  -TRT  ■  =  tAL  ■  I'* 

r  at  '■  -h. 

0  1  7 

ECU'  ISTRT  '=ECC ■  IS 

ThT  •  +  .  i.iCi ; 

0  1 75  730 

::=  istrt*.  14 

0 1 7? 

CRLL  PLOT'  VHT' 

ISTRT  -  .2  • 

0  1  U 

ISTR7= I STRT  + I 

0 ;  ! 

IF  'ISTRT  .LT.  lu 

i •  GOTO  710 

C  1  7  j 

IF  -'COUNT  .CT.  O' 

GOTO  700 

015  7  730 

CONTINUE 

0154 

CRLL  PLO 1 '0. ■ u.  • 

3  > 

i 


GO  HLIE 
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PROGRAM  DS2  (CONT) 


0 1  .::S 

WRITE' G. 17* 

01  c.1-’. 

LINE- 5 

013" 

DO  $00  1=1-1 00 

0 1 33 

LINE=LIME+l 

0  i  33 

IF  -LINE  .  LT.  5 >  001 0  750 

0 1  30 

WRITE' 0. 13  ■ 

0 1  3  i 

L I NE=0 

0132 

?50 

WRITE'  to  1$  ■  I.  TOT'  1«1  1  >  TtT', 

0 1 33 

REAR-  I  '--RBhR1  I  *  - R  hL  r 

013*; 

hTHETR'  1  'aflTHETrt'-  !  ■  120. 

0135 

BWGHT •'  I  ■  =BWGHT  •  I  '-LOCALE 

0133 

YD  AT  ft'  I  <  1  :'=BHGHT  •  I  1 

013? 

YDATRC I  -  L‘  *=R6hRi  I  * 

0 1 33 

YDRTR> I j  3>=ATHETA'  I  • 

0133 

$00 

CONTINUE 

0200 

DO  1000  I=li3 

020 1 

CALL  PLOT ■  .14.  YDn  r A'  1  •  I  • .  3  . 

0202 

DO  90O  I $TA T= 1 . 1 00 

0203 

:<=I$TAT*.  14 

0204 

CALL  PLOT'  YDATA'  I5TAT *  I  • 

02  05 

$00 

CONTINUE 

0203 

IF  I  .  EO.  1  ■  CALL  PLOT'  >..  . 

0207 

IF  ■'  I  .EO.  2'  CALL  F  LOP  u. . 

0203 

1  000 

CONTINUE 

0203 

CALL  URITE 

0210 

RETURN 

0211 

END 

0212 

END! 

J  •  S.'-ATHETA'  i 
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